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Abstract 

For any prime p the density of elements in F* having order, respectively 
index, congruent to a(mod d) is being considered. These densities on aver- 
age are determined, where the average is taken over all finite fields of prime 
order. Some connections between the two densties are established. It is 
also shown how to compute these densities with high numerical accuracy. 
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1 Introduction 

Let F* be the multiplicative group of a finite field Fg and let x G F*. The order of x 
is the smallest positive integer k such that x^ = 1, the index is the largest number 
t such that x^''-^)/* = 1. Note that t = [F* : (x)]. Let 7r(x) denote the number 
of primes p < x. Let 6{p; a, d) and p(p; a, d) denote the density of elements in 
F* having order, respectively index, congruent to a(mod d). It is not so difficult 
to show that both limx^oo N{a,d){x)/TT{x) and limx^oo N'{a,d){x)/'K{x) exist, 
where 

N{a, d) (x) = 5{p] a, d) and N'{a, d) (x) = p(p; a, d) 

p<x P^x 

These limits are denoted by 6{a, d), respectively p(a, d). In this note these quan- 
tities are investigated. The following result is characteristic of the type of results 
that will be established. (If a and b are natural numbers, then by {a,b), re- 
spectively [a,b], the greatest common divisor, respectively the lowest common 
multiple of a and b are denoted. By 7(a) = HpiaP square free kernel of a is 
denoted. Throughout the letter p will be used to indicate primes.) 

Theorem 1 

1) For every B > one has 

( X 

N{a, d){x) = 6{a, d)Li{x) + Ob 



log^ X 
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where 

oo oo / N 

S{a,d)= V V 77 

T^f mr(p[r['m,a\) 

(l + ra,cZ) — 1 (m,d)|a 

and t/ie implied constant depends at most on B. 
2) One has 

(p{{a,d)) ^ ^ ij,{m) 



(l + ra,d) = l {m,d) = l 

3) If d\di and 'y{di) ='y{d), then 6{a,di) = ■^6{a,d). 

4) One /ias 



5(0, t^) = ^ n ~~r and 2d) 

d , , 1 

5) For s > 1 one has 



p\d F 



|5(0, 2(i) z/ d is odd; 
5(0, 2(i) if d is even. 



r/ r)s\ _ / 2^ '^/S z/a zs even; 
^l«,2 j - l^i-s/a ^fa^sodd. 

6) Let q be a prime and q\ a. Then 

= 7 T^^r^ -T - ^— -p(--,g). 

(g — l)(g^ — 1) q^ — q — 1 a 

7) Pwt W'd(a) = {0 < r < c/ : (1 + ra, c/) = 1}. One has 

In the next subsection a characteristic zero version of 6{a, d) and p(a, d) will 
be discussed. Indeed, these characteristic zero quantities (exhibiting far more 
complicated behaviour) motivated the author to study 6{a,d) and p{a,d). The 
behaviour of these characteristic zero quantities turns out to have many resem- 
blances with that of 6{a, d) and p{a, d). 



1.1 Connections with characteristic zero 

Let g G Q\{ — 1, 0, 1} and p be a prime. By i^pig) the exponent of p in the canon- 
ical factorisation of g is denoted. If i^pig) = 0, then g can be considered as an 
element of F* with order ordg(p) and index rg{p). Let Ng{a,d) and Ng{a,d) de- 
note the set of primes p with i'p{g) = such that the order, respectively index of 
(/(mod p) is congruent to a(mod d). 

In case g = 2 the set N'g{a,d) was first considered by Pappalardi [1^], who 
proved that it has a natural density Pg{a, d) under the assumption of the Gener- 
alized Riemann Hypothesis (GRH). For the general case see [S]. 

The methods of Bj can be extended (see jTHI) to show that under GRH the 
set Ng{a, d) has a natural density 5g(a, d), the evaluation of which seems to be far 
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less easy than that of Pg{a, d). By 6g{d) the vector (5g(0, d), - ■ ■ , 6g{d — 1, d)) (if 
it exists) is denoted. Up to this century only Sg{2) had been evaluated. Recently 
Chinen and Murata [21 1121 computed ^£,(4) (on GRH) under the assumption that 
g is a positive integer that is not a pure power. In [H], on GRH, Sg{3) and 5g(4) are 
evaluated for each g e Q\{— 1,0, 1}. If d\2{a,d), then 6g{a,d) can be evaluated 
unconditionally, cf . [H [H IH 120] • 

Let G be the set of rational numbers g that cannot be written as —g^ or gQ 
with h > 1 an integer and go a rational number. By D[g) the discriminant of the 
number field Q(\/|^) is denoted. The functions N{a,d){x) and N'{a,d){x) can 
be considered as naive heuristic approximations of Ng{a, d){x) and Ng{a, d){x) (if 
S is any set of non-negative integers then S{x) denotes the number of elements in 
5* not exceeding x). For more on heuristics and primitive root theory, see [HEI- 
Theorem 121 and Theorem IHl say that as D{g) becomes large, the naive heuristic 
for 6g{a, d) and Pg(a, d) become more and more accurate. Some related numerical 
material is provided in Table 3 and Table 4. The next result is proved in pUj . 

Theorem 2 (GRH). If D{g) oo with g EG, then 6g{a,d) tends to 6{a,d). 

In this paper the following similar (but easier) result will be proved. 

Theorem 3 (GRH). If D{g) oo with g EG, then pg{a,d) tends to p{a,d). 

It turns out that both 6{a, d) and p{a, d) are much more accessible quantities 
than 6g{a, d), respectively Pg{a, d). In the light of the latter two theorems it thus 
seems of some importance to compute 5(a, d) and p(a, d), which is the purpose of 
this paper. The more complicated nature of Ng{a,d) versus Ng{a,d) is mirrored 
in the fact that S{a, d) is rather more difficult to compute than p{a, d). 

As a prelude to proving Theorem El its 'index equals t' analog is proved 
in Section |21 Theorem El is then proved in Section 3 (this involves evaluating 
p{a,d)). In Section 4 Theorem^ is considered. In Section 5 an Euler product 
Ay, involving a Dirichlet character x is studied and it is shown how p{a, d) and 
6 {a, d) can be expressed in terms of ^;^'s. Since can be evaluated with high 
numerical accuracy (Sections 6 and 7) this then allows us to evaluate p(a, d) and 
6 {a, d) with high numerical precision. 



2 Index t 

Let pp{t) denote the density of elements in F* having index t. Note that Pp{t) = 
f{{p — l)/t)/{p — 1) a p = l(mod t) and Pp(t) = otherwise. By the method of 
proof of TheoremjH (cf. (5| p. 161]) it is easy to show that, for every B > 0, 



p<x p<x ^ n=l ' \ / X o 

p = l{mod t) 

where Li(a;) denotes the logarithmic integral. Now 

^ tn^{tn) t^it) ^ n^{tn) " 
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The latter sum has as argument a multiphcative function in n. On applying 
Euler's identity, it is then inferred that 

ntipint) ' 

where 

A = T\{l- — i — - ) = 0.37395581361920228805472805434641641511 ■ ■ ■ 

p{.p-^)j 

denotes Artin's constant and 

r(t) = -T[ =^TT ^"^ 

Note that 

^ < K^) < -r^TTTT- (3) 



Combination of and yields the following result. 
Proposition 1 For every B > 0, one has 



^Pp(t) = Ar{t)U{x) + 



X 



B ^ ' 



log X 



Thus, the density of elements in ¥* with index t equals Ar(t) on average. 

Let Ng{t) denote the set of primes p with rg{p) = t. If g E Q\{ — 1, 0, 1}, then it 
can be shown [SlIIHl, under GRH, that Ng(t) has density 



n=l 



nt 1 



For an explicit evaluation of this density see [HI HHj. We can now prove the 
following result. 

Proposition 2 (GRH). Let g e G. If g > 0, set m = [2,D{g)]. If g < 0, 
set m = D{g)/2 if D{g) = 4(mod 8) and m = [4:,D{g)] otherwise. Put mi = 
m/{t,m). If g E G, then 

mg{t)) - Ar{t)\ < ^-^^ 



tmiLp(tmi] 



Corollary 1 (GRH). The density of Ng{t) exists and if D{g) — > oo with g EG, 
tends to the average density, Ar{t), of elements in F* having index t. 
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The latter corollary is the 'index equals f analog of Theorem |H1 



Proof of Proposition^ By and the evaluation of the deg ree [Q(Cfc, 9^^'') ■ Q] 
as given in [TH], it is deduced that 

m\kt m-[^\k 

On noting that ip{zw) > (p{w)ip{z), with w and z arbitrary integers and using 
that J2k < 2-21, the result then follows. □ 

Remark 1. The sum l/{k(p{k)) can be written as an Euler product of the form 
Hp -Pi (p) 7-^2 (p), with Fj{X) e Z[X] for j = 1,2 and monic. Using Theorem 2 of 
[H] such Euler products can be expressed in terms of values at integer points of 
the (partial) Riemann zeta-function. This enables one to evaluate these constants 
with hunderds of decimals of precision, see |T3|. A similar idea forms the basis 
of Theorem IHl and Theorem [3 



3 Computation of p(a, d) 

Equation (jH) suggests that, under GRH, one should have 
Proposition 3 (GRH). If g e Q\{-1,0, 1}, then 

fi{n) 



fea(mod d) t=a{mod d) n=l i^y^rit^d ^ ) • Q] 

(In this proposition and in the sequel sums over t are assumed to run over positive 
integers only.) Indeed by [12], cf. |H|, Proposition El is known to be true. Similarly 
one would expect that p(a, d) satisfies (|7j) as can indeed be proved. 

Theorem 4 For every B > one has 

N'{a,d){x) = p{a,d)U{x) +0 ( , (6) 

Vlog xj 

where 

p{a,d)=A y: r{t)= E ^n/^ (7) 

t=a(mod d) i=a(mod d) p\t 

and the implied constant depends at most on B. 
Proof. One has 

p<X t\p~l ^ 

t = a(mod d) 

On using that ip{n)/n = J2m\nf^(^)/^^ '^^^ obtains 

p{m) 



N'{a,d){x) = Y E E 



p<x t\p-l rn\E^ 

t = a(mod d) ' t 



mt 



Writing v = mt and bringing the summation over p to the inside, one obtains 



N'{a,d){x)= E ^ E -"(i) E 1- 

V<X — 1 t|" P<a; 

t = a(mod d) p^l(mod v) 

The summation range is spht up into the range v < log^"^^ x and log^"*"^ x < v < 
X. In the former range the Siegel-Walfisz theorem (see e.g. fr7\ Satz 4.8.3]) is 
invoked and for the latter range the trivial estimate ^p<j. p=i(mod 1 < is 
employed. Let d{v) denote the number of divisors of v. Together with the trivial 
estimate | "^^^ t=a(mod — '^('^) ^ which holds for every e > 0, one 

concludes (cf. jHl p. 161]) that © holds with 



p{a,d) = 



^-lt\v, t=a(mod d) 



Interchanging the order of summation and using ^ one infers that 

t=a(mod d) vi=l t=a{mod d) 

This concludes the proof. □ 

Let a > 0. As d becomes large, the first term in the second summation in ((7j), 
Ar{a), tends to be dominant by Corollary 0J In particular, lim^^oo ^(0, c?) = 
Ar{a). 

Proposition 4 One has 

j\ 1 J nj\ fp(0, 2d) if d is odd; 
(i9?((i) |^3p(0,2a) if d is even. 

Proof. Note that 

r((im) 



p(0,ci)=r(rf)5^ 



1 r{d) ■ 

Since r{dm)/r{d) is a multiplicative function in m, the identity for p{0,d) then 
follows on applying Euler's identity and noting that ip{d)/d = Y[p\d(^ ~ 1/p)- The 
identity for p(0, d) together with the observation that p(0, 2d)+p{d, 2d) = p(0, d), 
then yields the truth of the remainder of the assertion. □ 

By uj{m) the number of distinct prime divisors of m is denoted. 
Proposition 5 (GRH). Let g & G and m be as in Proposition^ then 



Pg{a, d) — p{a, d) 



< 



rrup[m) 
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Proof. By Theorem IHl Proposition 121 and (jSj) one infers on putting kt = v that, 
under GRH, 



Pg{a, d) — p{a, d) 



\p{k)\ 



- ^ ^ ktifikt) 

Tn\kt 



E 



i=a(mod d) IMi)l 



On noting that 

St|i), t=a(mod d) 1/^(^)1 ^ 

Vipiv) ~ ^-^ 

the result follows. 



_^ mvip{mv) ~ mip{m) ~ mip{m) 



□ 



Since 2^^"''^ / {rrnp{m)) tends to zero with increasing m and m tends to infinity as 
D{g) tends to infinity, Theorem El is a consequence of Proposition jSl 

The following result is concerned with Q-linear relations between the p(a, (i)'s 
with d fixed. 

Lemma 1 Let ai, ■ ■ ■ , a^(^d) be representatives of the reduced residue classes mod 
d. Then, for every integer a, 

p{a, d) G Q[p(ai, ci), ■ ■ ■ , p(a^(d)-i, d)]. 

Corollary 2 // d\di and Pi, - ■ ■ , are representatives of the reduced residue 

classes mod di, then 

Q[p(«i, c?), ■ ■ ■ , p{(^v{d)~i, d)] C Q[p(/?i, rfi), ■ ■ ■ , p(/?<^(rfi)_i, rfi)]. 
Proof of Lemma^ It is easy to show that 

Y^p{a,,d) = A Y: rit) = llll~----)eQ. (9) 

i=i (t,<i)=i pM 

It is thus enough to show that p(a, d) E Vd := Q[p(ai, d), ■ ■ ■ , p(«^(rf), rf)]. 

Let a = a/ {a, d) and S = d/ (a, c?). Note that (a, 5) = 1. Let 6i be the largest 
divisor of (a, d) with (5, 5i) = 1 and write (a, d) = S162. If a, 5, 61 and ^2 are 
being used for integers other then a and d, then this will be made explicit in the 
notation. Thus the meaning of Si{aj, d), which appears later in the proof, should 
be obvious. Note that 

p{a,d)=A Yl riSiS^t) = Ar{S2) ^ ^^'^i^)" 

t=ci(mod S) t=a(mod S) 

The proof proceeds with induction with respect to the number of distinct prime 
divisors of Si. If ^^(^i) = 0, then 61 = 1 and one has to show that p{a,6) G Vd, 
where a(mod 6) is a reduced residue class mod 6. Since a(mod 6) lifts to d/6 
reduced residue classes mod d, this is clear. If a;(5i) = 1, then 61 = with q a 
prime and e > 1. Then one has 

p(a,rf) = Ari62) Yl ^(^'^) 

t=a{mod (5) 
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= r{q')Ar{52) J] r{t) + ^r{52) J] r{qH) 

t=c(mod S) ^ t=^(raod 5) 

t^O(mod 9) 9 ' ' 

= Ci + -r(52) V r{qH) 

^ t=-|(mod <5) 

^ te^(mod 5) 

where Cj G V^. On choosing n > 1 to be such that = l(mod S), one infers that 
p{a, d) G Vrf. Suppose the result has been proved for all a and d with ^^(^i) < m 
for some m > 1. Then consider next a and ci with uj{5i) = m + 1. One has 

j = l t = ot{mo<l 5) 

t=j{mod 6-^) 

Note that 

A Y rm = A E r(52)r(5i(j,5i)t). 

t = cK(mod i5) t = Q(mod (5) 

t=j (mod ^i) f- . , fmod '^1 ) 

The latter sum equals a rational multiple times 

A Yl r{52)riij,5i)t) = p{aj,d), 

t = a(mod S) 

t=-7-rX r (mod . ^\ . ) 

for some integer Oj. Note that 8i{aj,d) = {j,Si). Thus by the induction hypoth- 
esis all terms in (fTUj) with ■y{Si) \ j are in Vd (where 7((5i) denotes the squarefree 
kernel of One infers that 

p(a,rf) = di + Ar{52) Y ^^^^^^ 

t=Q:(mod 5) 
t=0(mod 7(5i)) 

where 5i G V^. On choosing n > 1 to be such that 7(5i)" = l(mod 5) one infers 
that p(a, (i) G V^. □ 

Example 1. The result says that p(a, 6) G Q[p(l,6)]. Indeed, p(0, 6) = p(3, 6) = 
1/12. Furthermore p(2, 6) = 1/12 + 3p(l,6)/5, p(4, 6) = 1/3 - 3p(l,6)/5 and 
p(5,6) = 5/12-p(l,6). 
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4 Proof of Theorem U 



Proof of Theorem^ 1) Note that 

P<X r\p-l ^ 

^ =a{ mo d d) 

Proceeding as in the proof of Theorem |5 one infers that 

rKx—1 fji-^ x — l P<x, p=l(mod rm) 

^ p = l + ra(mod rd) 

Now for the inner sum to be non-zero the two congruences must be compat- 
ible. By the Chinese remainder theorem this is the case if and only if 1 = 
1 + ra(mod r{m,d)), that is if and only if a = O(mod {m,d)). If the two con- 
gruences are compatible, then they form a reduced residue class if and only if 
(1 + ra, d) = 1. If the residue class is not reduced it contains at most one prime 
and the contribution of these primes to N{a,d){x) is bounded in absolute value 
by E„<x d{v)/v = O(log^x). 

The summation range mr < a; — 1 in (fTT|) is split up into the range r[m, d] < 
log'"x and the range r[m,d] > log*" a;, where C is to be chosen later. All error 
terms arising in this way are easily seen to be of the claimed order of growth, 
except the error term 

r mr,f(r[m,d\) 
which arises on completing the sum 

|/i(m)| 



Li(-) E E 



{l + ra,d) — l ^ ^ 

r[m, d]<log^ X 



mr(p{r[m, d]) ' 



to S{a, d). On noting that r[m, d\ > log*" x implies rmd > log*" x and using that 
ip{zw) > ip{z)ip{w), one obtains, cf. the proof of part 2, 



E(x) = O 



^ Li(x) V - |/i(mi)| |/i(w22)| ^ 

\ rm2>log^ x/{dm^) / 



From this E(x) is easily seen to be Ob{x/ \og^ x), when C is chosen to be suffi- 
ciently large. 

2) By part 1 it is enough to show that 

/i(m) (p{{a,d)) /i(m) 



{m,d)\a (m,d) — 1 
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Note that 



m.i\{a,d) ™=l T-v mi ^ 



On writing m = mim2 one obtains 

/i(mim2) 



rm\{a,d) ^ "2=1 i ^ v-v ^ ; 

>^ /x(mi) fJ'{m2) 
^ m\ ^ rm2^(m2rd) 

mi\(a,d) "12 = 1. (m2,d/mi) = l 

{™2.'"l) = l 



(m2,d) = l 

V9((a,d)) ^ /i(m) 



l*^) ^^rf rm(p(mrd) 

(m,d) = l 

3) The condition on d and ensures that (1 + ra,di) = 1 iff (1 + ra,d) 
1 and {m,di) = 1 iff {m,d) = 1. Furthermore one has ip{{a,di))/{a,di) 
(p{{a,d))/{a,d). By part 2 one then finds 

(p{{a,d)) ^ ^ /i(m) 



(a,c?) ^ frt mnp{mrdi)' 

{l + ra,d) = l {m,d) = l 

On noting that ip{mrdi) = (p{mrd)di/d, the proof of part 2 is completed. 
4) By part 2 one has, on writing mr = v, 



Am - Emit,, {m,d)=l 



d ^ v(p(vd) 

v=l ^ ' 

The inner sum equals one if 7(f) |(i and zero otherwise. Thus 

'l{v)\d ' l(y)\d p\d P 

The formula for 6{d, 2d) easily follows from that of 5(0, d) and the observation 
that (5(0, 2d) + 6{d, 2d) = 5(0, d). 

5) An easy consequence of part 3 and part 4. 

6) Using part 2 and Q, one infers that 

5(a,g) = 5(0, g)- f ^^^^ 



— ^ mt(j)(mtq) 

t=_i(mod g) ™=l 

a ^ (m,g) — 1 



E E 



fi{m) 



q — 1 ^ ^ ■mtMmt) 

t=-l(mod g) ^ ^ 

a V ^) (Tn,g) — 1 
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te-l(mod q) 9(9-1) 

The proof is then completed on invoking part 3 and ((7|). 
7) By part 2 one can write 

5(a d) = '^^^^'^^^ 1 /i(m)y(rrf) 

' (a, d) ripird) mLp{mrd) 

(l + ra,d) = l (m,d) = l 

Denote the inner sum in (jl2p by One has 

= n(i-^)n(i-^) 

p\d 

A 1—1 1 — 1 

^ ^ TT TT p(p-i) 

n (1 ^ ) * M ^ J-J- 1-^ ■ 

llp|dV p(p-l)'' p|r. P{P-1) p\{r,d) 

One thus obtains that 



1 



(<■, -i) n,i,(i - .,t:w ^. 1 - ? 



lij — ci(mod d) 

Note that 



This equation together with the latter one derived for 5(a, d) and 0, then yields 
the result. □ 



5 The densities and 

Given a Dirichlet character x mod (i, let = x * that is /i^ denotes the 
Dirichlet convolution of x and the Mobius function /i. Note the following trivial 
result. 

Lemma 2 The function hy. is multiplicative and satisfies h^{l) = 1 and further- 
more h^lp"^) = x{pY^'^[x{p) " l]j where the convention that = 1 zs adopted. 

In particular if x is the trivial character mod d, then 

I- otherwise. 

By using one of the orthogonality relations for Dirichlet characters, the following 
result is easily inferred. 
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Lemma 3 Let a(mod d) be a reduced residue class mod d. One has 

t\v 

where x runs over the Dirichlet characters modulo d. 

The reader is referred to Section 2.4 of ^8J for some further properties of h^. 
In what follows sums of the form 

{v,d)=l ' 

will feature. It is easy to see that this sum is absolutely convergent. Since its 
argument is multiplicative, one then obtains that the latter sum equals 

Note that if x is the principal character, then = 1. For a fixed prime p and 
a G M, < a< 1 let 



(p2 _ e27riQ)(p_ 2 

A tedious analysis shows that |/p(«)| as a function of a is decreasing for < a < 
1/2 and increasing for 1/2 < a < 1. Thus 

where the lower bound holds true iff a = 1/2 and the upper bound iff a = 0. It 
follows that \A^\ < 1 with = 1 iff x is the principal character mod d. 

If x' is the primitive Dirichlet character associated with x, then the Euler 
products of A^ and A^/ differ in at most finitely many primes and hence can be 
simply related. 

It will be shown in Theorem El that p{a,d) and 6{a,d) can be expressed in 
terms of A'^s, where x ranges over the Dirichlet characters mod d. The proof 
makes use of the following proposition. 

Proposition 6 Let a > 1. One has 



2 _|_ x(p)-i 



p^-x(p) 



_l_ (x(p)-i)p 

X(mod <5) Pl» (P^-X(P))(P-1) 



where w = {a,d), a = a/w and 6 = d/w. In particular, if (a, d) = 1 then 
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Corollary 3 If '-f{{a,d)){a,d)\d, then 

wipiw 



p(^a d) = P^^'^'^ TT _ ^ ~ ^ 



p\{s,w) p(p-i) 

Proof of Proposition |B] From (jH)) and Lemma El one easily infers that 
pia d) — ~)J7)!)'S^ ^^^'"^'^^'^^ 

' ' x(mod 6) v=l > 

On noting that the argument of the inner sum is multiplicative in f , the result 
follows on applying (fT^ and Euler's product identity. □ 

Example 3. One has p(0, rf) = l/{dif{d)) (in agreement with Proposition HJ. Let 
Xs and Xi denote the non-trivial character mod 3, respectively mod 4. One finds 
p(±l,3) = 5(1 ±A^3)/12 and p(±2,8) = 5(1 ± A^J/12. Let x be the character 
mod 5 uniquely determined by x(2) = i- One has 



19 



p(3,5) = -(l + 2Re(2AJ-A, 



.2 , 



80 ■ - ^ 

Using Table 2, these densities can then be numerically approximated. 

Example 4. One has 

p(2,8) = 3p(l,4)/4 rp(3,9) = 8p(l,3)/45 / p(2, 12) = 3p(l, 6)/4 

p(6,8) = 3p(3,4)/4 ^'"'^ \p(6,9) = 8p(2,3)/45 ^""'^ \ p(10, 12) = 3p(5, 6)/4. 

One can now infer how the densities can be expressed in terms of Ay_^s. 
Theorem 5 Let a and 5 he as in Proposition IHl Then 

p{a,d) e Q(Cord„(5))[^xir--,^X^(,)], 

where Xi? ' ' ' ? XipiS) (^"i^^ the characters mod 6. 

Let A denote Carmichael's function, that is X{d) equals the exponent of the 
group (Z/dZ)*, then 

5(a,rf) GQ(CA{d))[Axir--,^x.(.)]- 

Proof. The first part is a straightforward consequence of Proposition ^ The 
second part follows on applying part 7 of Theorem Q and Proposition [HI together 
with the observation that if 6\d, then any character x' mod 6 can be lifted to a 
character x mod d, such that A^' = cA^, where c G Q{(x{d))- ^ 

The next result follows on combining Proposition IHl with part 7 of Theorem ^ 
Proposition 7 Suppose that q is a prime and q\ a. Then 
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The Euler product can also be expressed in terms of p{a, dys. This yields 



Thus, by ^ and p(a, d) > 0, one finds 



with equality iff x is the principal character mod d. 

6 The numerical evaluation of 6(a,d) and p(a,d) 

Consider the numerical evaluation of the constant A^. To this end it turns out 
to be more convenient to consider 



p V [p^-x(p)](p- 1)7 -Ji^ V p(p-i) 

Recall that L{s, x), the Dirichlet series for the character x, is defined, for Re(s) > 
1 by L{s,x) = T.n=iX{n)/n'. 

Theorem 6 Let pi{= 2),p2, - ■ ■ be the sequence of consecutive primes. Let x be 
any Dirichlet character and n > 31 {hence pn > 127). Then 



n , 

= R^AL{2, x)m XW, X)\{[1 + 

k=l ^ 



x{p) \ /I x{Pk).,^ x{Pk) 



)(1 



Pkipl - Pk - 1) J pI pI 



with 

1 , , 1 

< \Ri\ < 1 + 



+ Pn+1 Pn+1 

Proof. The first step is to note that 



oo ^ 

B^ = AL{2,x)H3,xW,x)l[U 

k=l ^ 



x{p) \ n x{pk).,^ x{p, 



(1-^)(1 



•k 



Pk{pl - Pk - 1) J pI pI 

An upper bound for the kih term in the latter product is given by 



5 (2 + 2t + + t^) 
1-t-t^ ' 



where t = l/pk- For t > 127 some analysis shows that the latter expression is 
bounded above by 1 + 1^-*^. Using this one obtains 

-p-r / 1 \ 1 I dt 1 



P>Pn ^ ' m>p„ 
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A similar argument allows one to obtain the lower bound. 



□ 



Since the Artin constant (see e.g. and L(2, x), L(3,x) and L(4, x) can be 

each evaluated with high numerical accuracy, Theorem IHl allows one to compute 

with high numerical accuracy. Using Proposition IHl and part 7 of Theorem Q 
p(a, d), respectively S{a, d), can then be evaluated with high numerical precision. 

A more straightforward, but numerically much less powerful, approach in 
computing p{a,d) and 6{a,d), is to invoke part 7 of Theorem ^ and compute 
p{a,d) using the identity p{a,d) = AJ2t=a{mod d)''^i^)- One has the following 
estimates. 



Proposition 8 Let x > Q. One has 



< p{a,d) - A ^ r{t) < 



1.28 



t — a(mod d) 
t<x 



X 



Corollary 4 Let a > and a + d > Q, then 

1.28 



< p(a, d) - Ar{a) < 



a + d 



The most important ingredient of the proof will is the following lemma (the idea 
of which was suggested to the author by Carl Pomerance [IE])- 



Lemma 4 For x > 6 one has 



^ 1 2.1 

nip{n) X 



n>x 



Proof. Using that f{n) > log(2n)/(?2log2) and that ^fc>yl/A:^ < 1.075/?/ for 
y > Q, one finds, for x > 6, 

\p{d)\ 



^ tM 5Z ^2 5Z ^(d) 

t>x ^ t>x d\t ' 



^ d^ioid) ^ 

d=l ' r>x/d 

^ 2^ d^o,(d) X ^ d^id) 

^ ^ rf31og2 X ^ dip(d) 
d>x/6 ^ d<x/6 ' 

< C(2) r i^^, + H[^C(2)C(3) 
- ^^'V[./6]tMog2^'+ X C(6) 



< 



< 



< 



C(2) (21og(2[x/6]) + l) 1.075 C(2)C(3) 
4 [x/6]'log2 X C(6) 
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For X > 45000 the latter upper bound is bounded above by 2.1/ x. After calcu- 
lating l/l'^V^l'^)) with enough precision (see Remark 1) and using that 



t^^^H tt^^(^) ^i^^^H ' 

the result follows after verification in the range (6, 45000) (this verification is 
easily seen to require only a finite amount of computation, cf. [3 Lemma 4]). □ 



The above argument can be easily adapted to show that 

1.943 

^ n^{n) C(6) 



n>x 

Proof of Proposition^ From ((7j) and r{t) > 0, one infers that 

A Yl r{t)<p{a,d)< Yl r{t) + Aj2rit) 



By one has 



t — a(mod d) £ — a(mod d) t^X 



t>x t>x ' 



On invoking Lemma El the result then follows. □ 

The terms AL{2, x)L{3, x)L{4, x) in Theorem |H1 form the beginning of an expan- 
sion of in terms of special values of L-series. 

Theorem 7 Define numbers G^^}_^ by 



r 2-^ (^1 _ _ ^2dy/d 2-^ i+1 



d\r ^' " " ' J=0 



One has 



^ ^ ^ ) r=l fc=3r+l 

where {-iy-^\{k,r) = 0^^,+^ G Z>o. 

Note that as formal series {1 — z — z"^)^^ = J^JLo-^j+i^'' where Fj denotes the 
jth Fibonacci number (thus Fq = 0, Fi = 1 etc.). The numbers defined by 
{1 — z — z^)^''' = X]jlo -^j+i^"' known as convolved Fibonacci numbers and 
hence a reasonable term for the integers G^+i might be 'convoluted convolved 
Fibonacci numbers'. For the convenience of the reader Table 1 gives a small 
sample of these numbers. 

The positivity of the numbers (— l)^'^^A(/c, r) is established in [U], where the 

(r) 

numbers Gj^^ are investigated. The argument uses Witt's dimension formula for 
free Lie algebras. The remaining part of Theorem [3 follows from the following 
more general result. 



oo oo 
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Theorem 8 Suppose that f{z) allows a formal power series in z having only inte- 
ger coefficients, i.e. f{z) = Ej>ia(j>^ with a{j) G Z. Let g{z) = Y.j>i 
and let jo > denote the smallest integer such that a{j) 7^ 0. Let 



00 

d\r j=0 

Then, as formal power series in y and z, one has 



r)z^ . 



00 00 



i-yf{z) = l[ n (1 - ^V)'^^''\ (14) 

k=lj=kjo 

Moreover, the numbers h{j, k) are integers. 

Let e > be fixed. The identity fll4|) holds for all complex numbers y and z 
with g{\z\)y < 1 — e and \z\ < p^, where pc is the radius of convergence of the 
Taylor series of g around z = 0. If, moreover, > 1/2, g{l/2) < 1 and X]p5'(p) 
converges, then 



00 00 

hij,k) 



i-x(p)/(-))=nn^(^'^ 

p k=l j=kjo 

Proof of Theorem^ Note that 



(15) 



:i - Yx' 



l-J^ I V l-X-X 



i-x 



By the first part of Theorem |S1 one infers that, as formal series, 

V V3 00 00 

(1 + , _y_ = (1 - yxr\i - Y'x') n n (i-^^n^-^^^^ 

r=l fc=3r+l 



fc-3r+l 



on noting that G^p = 1 for r = 2 and G^'^^ = for r > 2. Apply the second 
part with f{z) = —z^/{l — z — z^) (and hence g{z) = z^/{l — z — z^^). The 
Taylor series for g has radius of convergence pc = (v^- l)/2 > 0.5. Note that 
g{x) < 1/2 for all < x < 1/2. Furthermore, Ep^(l/p) < Ep |i < 00. □ 

Remark 2. The Dirichlet character is an example of a completely multiplicative 
function h, i.e. h{nm) = h{n)h{m) for all natural numbers n and m. If one 
defines L{s, h) by L{s, h) = Yl'^=i h{n)n~^, then under the same conditions, one 
may replace x i^i Theorem |H1 by any completely multiplicative function h satisfy- 
ing \h{n)\ < 1. 

Remark 3. If < pc < 1/2 and ^p^i/p^gil/p) converges, then an identity of 
the type (fT3j) still holds, but with Dirichlet L-functions being replaced by partial 
Dirichlet L-functions. The idea is just to leave out the local factor 1 — x(p)/(l/p) 
for sufficiently many small primes p and then proceed as before, cf. the proof of 
Theorem 1 of [Hj (in the formulation of Theorem 1 there, replace pno + 1 > 1//5 
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(a typo) by p^o+i > P). 



Remark 4. The conditions in the latter part of the theorem ensure that f{z) = 
O(z^) for small z. This ensures on its turn that in the double product in |T3|) 
only factors L{j, x^) with k > 1 and j > 2k > 2 appear. 

The proof given here of Theorem [7| rests on the following lemma. 

Lemma 5 Suppose that f{X,Y) = Ylj k'^U^ ^)-^''^^ with a{j,k) integers and 
/(0,0) = 0. Then there are unique integers e{j,k) such that, as formal series, 
one has 

oo oo 

1 + f{X, Y) = Y[ Yl (1 - X^Y'')<^'''\ 

)=0 k=0 

(j,fc)#(0,0) 

Proof. The term X^^Y^^ is said to be of lower weight than X^^Y^"^ if ki < k2 or 
ki = k2 and ji < j2- Suppose that X^Y^ is the term of lowest weight appearing 
in f{X,Y). Then consider (1 + f{X,Y)){l - X^Y'')-''^^'''\ This can be written 
as 1 + g{X, Y) where all the coefficients of g{X, Y) are integers and the term of 
lowest weight in g{X, Y) has strictly larger weight than the term of lowest weight 
in f{X,Y). Now iterate. 

It is not obvious from this argument that if one starts with a different weight 
ordering of the terms X^Y'^ we end up with the same integers e{j,k). Sup- 
pose that h{X) has integer coefficients, then the coefficients e(n) in 1 + h{X) = 
n^i(l ~ X^Y^"-^ are unique, cf. 0. Hence, by setting X = 0, respectively 
Y = 0, one obtains that e{0,k), respectively e(j, 0) are uniquely determined. 
Setting F = X™ one obtains that 1 + /(X, X™) = Un=ii^ ' X")^("), where v{n) 
is uniquely determined and v{2m) = e(2m, 0) + e(m, 1) + e(0,2). The unique- 
ness of e(0,2), e(2m, 0) and /(2m) then implies the uniqueness of e(m, 1). The 
proof will be completed by using induction. So suppose one has established 
that e(j, k) with k < r for some r > 1 are uniquely determined. Using that 
f((r + 2)m) = ^fc=o^(('" + 2 — k)m, k), one infers by the induction hypothesis 
and using that e(0, r+2) and v{{r+2)m) are uniquely determined, that e(m, r + 1) 
is uniquely determined. □ 

Proof of Theorem\^ By Mobius inversion and the definition of H^'^\z) one infers 
that 



oo 

d\r d\r j=0 



from which it is inferred that 



oo oo 

kd 



r=l k=l j=0 d=l 

The latter identity with both sides divided out by y can be rewritten as 

f{z) ^^h{j,k)kz^y^~^ 



sr^ sr^ ityj, I 
l-yf(z) 2^2^ l-z^y^ 
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Formal integration of both sides with respect to y gives 



oo oo 



- log(l - yf{z)) = -Y.Yl ^(^■' ^) - ^'^/')' 

k=l j=0 

whence 

oo oo 
k=l j=0 

On writing fi{z) = f{z)/z^° and yi = yz^" and expanding l — yifi(z){= l—yf{z)) 
in terms of yi and z, it is then seen that (fT^ holds. The integrality of h{j, k) 
follows by Lemma El 

The formal argument can be certainly made rigorous in the situation where 



oo oo 



5^X;iMj,fc)A:^Vl<oo, (16) 

^==1 3=30 

where one is in the situation of absolute convergence and interchanges in order 
of summation are hence allowed. Note that for x > 0, (7 is an non-decreasing 
function of the real variable x. Now note that 



00 

Y,\Kj.r)z^\<-Y.3i\A'r"<9{\Ay 



r 

j=jo d\r 

where the assumption that \z\ < Pcis being used. The double sum in flTH|) is thus 
majorized by '^'^i'rg{\z\Y\y\^ which in the given {y,z) region converges. 

By a similar argument the convergence YlpYl'^^i'^^di^/pY ^ which is a conse- 
quence of the convergence of J2p9i^/p) (^^^ ^^^^ heie that g is non-decreasing 
as a function of the real variable x for x > and that g{l/2) < 1), ensures the 
convergence of the triple product 

n(i-x(p)/(i))=nnn(i-^)'*"- 

P ^ ^ ^ P k = lj=kjo ^ ^ ^ 

(From the theory of infinite products use that a product ^(1 + e^) is called ab- 
solutely convergent if ^ is absolutely convergent and that in an absolutely 
convergent product the factors can be reordered without changing its value.) On 
bringing the outer product over the primes p to the inside and using the Euler 
product for a Dirichlet L-series, the result then follows. □ 
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7 Tables 



Explanation to Table 1. Table 1 gives some values of convoluted convolved Fi- 
bonacci numbers G^^. These numbers are defined in Theorem [7| 

Explanation to Table 2. For every character x of modulus < 12, A^^ can be de- 
duced from the table below. In every case the value of x is given (in at most 
two arguments) such that x is uniquely determined by this. If x itself is not in 
the table, its complex conjugate x will be (in which case one has = A^) or 
X is the principal character (in which case A^ = 1). Although A^ for x not a 
primitive character can be easily related to A^/ with x' ^ primitive character, 
for the convenience of the reader the numerical approximations to A^ for the 
non-primitive characters are listed as well. 

Explanation to Table 3. An entry in a column having as header the number a 
and in a row starting with an integer d, respectively a — , gives the first five 
decimal digits of 6{a, d), respectively 6{a + 6, d). If an entry is in a row labelled 
~, let 6{a, d) be the entry directly above it. Then the number given equals 
N_ig{a,d){x)/Tr{x) with x = 2038074743 (and hence 7r(a;) = 10^). 

Explanation to Table 4. Similar to that of Table 2 (and with the same value of 
x). In case d = oo one has 6{a,d) = Ar{a) for a > 1 and Ng{a, oo){x) denotes 
the number of primes p < x with Vp{g) = such that g has index equal to a. 
Here x = 1299709 (and hence 7r(1299709) = 10^). 

Acknowledgement. The author thanks Yves Gallot for writing a C'^'^ program 
that was used to create Tables 3 and 4. Furthermore, he thanks the referee for 
some helpful comments. 
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Table 1: Convoluted convolved Fibonacci numbers Gj^ 



r\j 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


1 


1 


1 


2 


3 


5 


8 


13 


21 


34 


55 


89 


2 


1 


1 


3 


5 


11 


19 


37 


65 


120 


210 


376 


3 





1 


3 


7 


17 


37 


77 


158 


314 


611 


1174 


4 





1 


3 


10 


25 


64 


146 


331 


710 


1505 


3091 


5 





1 


4 


13 


38 


102 


259 


626 


1457 


3287 


7224 



Table 2: Numerical evaluation of A-^ 



d X 


/V 




3 


X(2) = -1 




-|- u. i ( oy ( ( izz4izyDo4i: ■ ■ ■ 


4 


X(3) = -1 




-T U.D^ODOUD / yODZOZO • • • 


5 


X(2) = t 




-r U.OD'iDoyDZD^l: / oOol • • • 








-l-i'n 00/10/11 OQ/I /1 0/1 7'?S 


5 


X(2) = -1 




-|- u.izyou / yooozoUoU • • • 


6 


X(5) = -1 




_L n 8(^Q88c:(^1 01 /I 81 71 

-r u.ouyooouiZi^toi < i ■ ■ ■ 


7 


X(3) = e''*'-^ 




_L n 01 87fiQ0Q8/19Q'^fiQ 

-|- u.zio ( oyzyo^tzyooy ■ ■ ■ 


' 






-LiD O'iE^/l 1 8/1 c;fifi7Q 

-rtu.zoo^tio^oooooo ( y ■ • • 


7 


X(3) = e*''*'-^ 




_i_ n 01 0R1 07fin/i7c;nR0 

-|- U.ZiZDiZ/ 01)4 / OUDZ • • • 








iTl 1 /1 1^1 SSQSRQflfiRI n 
— ZU.i40iooyoDyUoDiU • • • 


7 


X(3) = -1 




-1- n fil 1 '^oAA'^OQi (yvi'x . . . 
u.u i ioz440zy iyo io • • • 


8 


X(3) = 1 


Xl^j — -J- 


-1- n 8'?7QQ8p;n'^i oo'^fin . . . 
-p u.oo ( yyoouoizyoou • • ■ 


8 


X(3) = -1 


Xl^J — 


-p u.u^oijouij ( yoozozo • ■ ■ 


8 


X(3) = -1 


X(,oj - -i 


-J- n fin'?Qn78p;fi9fi7i fi7 . . . 
u.uuoyu ( oouzu ( ±u i • • • 


9 


X(z) = e 




_i_ n 1^7881 c:01 1 fiQOOO/l 

-i- u.o (ooioyiiDozyz4 • • • 








^(■(1 x\ 1 iriR i inn inoo^i . . . 

\l\j.OO 1 lUOl V\J\J L\)LoO 


9 


X(2) = e^'^*/^ 




+ 0.250710892521489- ■■ 








-i0.207858981269346- ■ ■ 


9 


X(2) = -1 




+ 0.173977122429634- ■ ■ 


10 


X(3)= I 




+ 0.779414790379699 ■ ■ ■ 








+iO. 123970019579663- •• 


10 


X(3) = -1 




+ 0.646539692640401 • • • 


11 


X(2) = e'^^/^ 




+ 0.657644343795360 ■ ■ ■ 








+i0.151998116640767--- 


11 


X(2) = e^^^^/^ 




+ 0.373259555803500 ■ ■ ■ 








+i0.208638808901506--- 


11 


X(2) = e^^^'/^ 




+ 0.187051722258759 • • • 








+i0.232381723173172 • • • 


11 


X(2) = -1 




+ 0.184204262987186- ■■ 


12 


X(5)= 1 


x(r) = -1 


+ 0.919500970946465- ■■ 


12 


X(5) = -1 


X(7)= 1 


+ 0.869885612148171 • • • 


12 


X(5) = -1 


X(7) = -1 


+ 0.841259078358102 • • • 
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Table 3: S{a, d) and approximation to (5_i9(a, d) 



(I 





1 


2 


3 


1 


5 


d^2 


0.66666 


0.33333 












0.66667 


0.33333 










3 


0.37500 


0.35599 


0.26900 










0.37502 


0.35602 


0.26897 








4 


0.33333 


0.16666 


0.33333 


0.16666 








0.33334 


0.16664 


0.33333 


0.16669 






5 


0.20833 


0.23542 


0.17799 


0.23400 


0.14424 






0.20831 


0.23572 


0.17829 


0.23373 


0.14395 




6 


0.25000 


0.06067 


0.12134 


0.12500 


0.29532 


0.14766 




0.25001 


0.06067 


0.12132 


0.12501 


0.29534 


0.14765 


7 


0.14583 


0.15968 


0.15483 


0.11905 


0.16351 


0.15567 




0.11581 


0.15965 


0.15167 


0.11915 


0.16367 


0.15573 




0.10141 














0.10129 












8 


0.16666 


0.08333 


0.16666 


0.08333 


0.16666 


0.08333 




0.16667 


0.08332 


0.16664 


0.08335 


0.16667 


0.08332 




0.16666 


0.08333 












0.16669 


0.08334 










9 


0.12500 


0.11866 


0.08966 


0.12500 


0.11866 


0.08966 




0.12501 


0.11866 


0.08966 


0.12501 


0.11868 


0.08965 




0.12500 


0.11866 


0.08966 










0.12500 


0.11867 


0.08965 








10 


0.13888 


0.07196 


0.14393 


0.08172 


0.06810 


0.06944 




0.13888 


0.07197 


0.14408 


0.08159 


0.06783 


0.06944 




0.16345 


0.03405 


0.15227 


0.07613 








0.16374 


0.03421 


0.15214 


0.07612 






11 


0.09166 


0.09890 


0.09811 


0.09904 


0.09848 


0.07170 




0.09166 


0.09889 


0.09805 


0.09904 


0.09859 


0.07180 




0.09910 


0.09303 


0.09297 


0.09523 


0.06113 






0.09939 


0.09303 


0.09297 


0.09526 


0.06133 




12 


0.12500 


0.03033 


0.06067 


0.06250 


0.14766 


0.07383 




0.12500 


0.03033 


0.06065 


0.06251 


0.14767 


0.07382 




0.12500 


0.03033 


0.06067 


0.06250 


0.14766 


0.07383 




0.12501 


0.03035 


0.06067 


0.06249 


0.14767 


0.07383 
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Table 4: p{a,d) and approximation to ^65537(0, (i) 








1 


2 


3 


1 


5 




0.50000 


0.50000 












0.49994 


0.50006 










3 


0.16666 


0.48915 


0.34417 










0.16662 


0.48924 


0.34414 








4 


0.12500 


0.41091 


0.37500 


0.08908 








0.12497 


0.41097 


0.37497 


0.08909 






5 


0.05000 


0.44143 


0.31320 


0.10036 


0.09498 






0.05000 


0.44150 


0.31322 


0.10035 


0.09494 




6 


0.08333 


0.38955 


0.31706 


0.08333 


0.09959 


0.02710 




0.08330 


0.38966 


0.31705 


0.08331 


0.09958 


0.02709 


7 


0.02380 


0.40253 


0.29923 


0.08966 


0.08471 


0.03881 




0.02380 


0.10263 


0.29923 


0.08962 


0.08170 


0.03881 




0.06123 














0.06122 












8 


0.03125 


0.38569 


0.30818 


0.07380 


0.09375 


0.02521 




0.03124 


0.38577 


0.30818 


0.07380 


0.09372 


0.02521 




0.06681 


0.01528 












0.06679 


0.01529 










9 


0.01851 


0.39347 


0.29075 


0.08696 


0.07829 


0.02983 




0.01851 


0.39356 


0.29075 


0.08694 


0.07829 


0.02981 




0.06118 


0.01738 


0.02358 










0.06117 


0.01740 


0.02358 








10 


0.02500 


0.38063 


0.30067 


0.07141 


0.08456 


0.02500 




0.02500 


0.38071 


0.30068 


0.07140 


0.08452 


0.02500 




0.06080 


0.01253 


0.02895 


0.01041 








0.06079 


0.01254 


0.02895 


0.01041 






11 


0.00909 


0.39040 


0.28866 


0.07722 


0.07791 


0.02698 




0.00910 


0.39047 


0.28865 


0.07721 


0.07791 


0.02698 




0.05513 


0.01768 


0.02331 


0.01118 


0.01909 






0.05541 


0.01769 


0.02331 


0.01419 


0.01908 




12 


0.02083 


0.37819 


0.29216 


0.07231 


0.07926 


0.02170 




0.02080 


0.37827 


0.29215 


0.07230 


0.07926 


0.02169 




0.06250 


0.01136 


0.02489 


0.01101 


0.02033 


0.00540 




0.06250 


0.01139 


0.02491 


0.01101 


0.02033 


0.00539 
















00 


0.00000 


0.37395 


0.28046 


0.06648 


0.07011 


0.01889 




0.00000 


0.37367 


0.28124 


0.06646 


0.06913 


0.01885 




0.04986 


0.00893 


0.01752 


0.00738 


0.01417 


0.00340 




0.04962 


0.00915 


0.01796 


0.00745 


0.01449 


0.00359 
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